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Using the BF version of pure Yang-Mills, it is possible to find a covariant rep- 
resentation of the 't Hooft magnetic flux operator. In this framework, 't Hooft's 
pioneering work on confinement finds an explicit realization in the continuum. Em- 
ploying the Abelian projection gauge we compute the expectation value of the mag- 
netic variable and find the expected perimeter law. We also check the area law 
behaviour for the Wilson loop average and compute the string tension which turns 
out to be of the right order of magnitude. 
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1 Introduction 



The problem of quark confinement in continuum QCD is still unsolved. Confinement of 
electric charges is observed if one is willing to use the lattice as a regulator but direct 
computations in the continuum theory are not easily performed due to the large scale 
nature of the problem at hand. 

Among the different explanations put forward in literature, we find that the so called 
dual Meissner effect is the most appealing and complete: if the QCD vacuum behaves as 
a magnetic superconductor color electric flux lines will form long and narrow tubes with 
color electric charges at its ends. The story of this idea is quite long and has received 
contributions from many authors j^, ^ . In this paper we will stick to the version of 
't Hooft in which all these ideas are brought to a high degree of deflnitiveness P, P]. 

More in detail, the whole idea of the dual superconductor can be phrased in terms 
of 't Hooft order-disorder variables. These variables are non-local quantum operators, 
W{C), M{C'), creating thin color electric and magnetic fluxes, along the curves C,C', 
in the QCD vacuum. In their turn, the curves C, C border two surfaces S, S'. In the 
absence of massless particles, these surfaces may become dynamical and we may have 
linkage when one of these curves (say C for example) intersects the surface (say S) 
spanned by the other curve (C). The linkage costs energy per unity of area and, roughly 
speaking, the associated effective action goes as the area. A, of the dynamical surface: 
in our previous example it will be < W{C) >~ ^'^'^^^■M^) ^ p^\\ these properties are 
summarized by the commutation relation of W{C), M{C') which is of braiding type and 
which contains the entire physics of the problem. Which surface is intersecting which 
curve is not known in this picture which supports all possibilities: Higgs phase. Coulomb 
phase, partial Higgs with conflnement and conflnement phase ||^. To support the picture 
of the dual superconductor the QCD vacuum must be flUed by magnetic vortex lines with 
non-trivial topology, generated by means of a true or effective dual Higgs dynamics. 

Despite the appeal of this quahtative picture one may still be left with a feeling of 
uneasiness since quantitative computation in the continuum theory have not yet been 
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performed. 

In order to expose the topological content of the theory, one of us has proposed to 
study pure Yang-Mills theory using the first order formalism (BF-YM). The main advan- 
tage of this formulation is that a color magnetic operator can be easily guessed in the 
continuum . In this formalism we show that our magnetic operator is the dual 't Hooft 
observable and compute its vev, employing the Abelian projection gauge |], finding 
the expected perimeter law. 

The first order form of pure Yang-Mills is described by the action functional 

2 

Sbf-ym = J Tr(z5 AF + ^BA*B) (1.1) 

where 



F = ]^F^^T^dx^ A dx" = dA + i[A, A] , (1.2) 



D = d + i[A, ■] and B is a Lie valued 2-form. The generators of the SU{N) Lie algebra 
in the fundamental representation are normalized as TrT"T^ = 5afe/2 and the * product 
(Hodge duality) for a p form in d dimension is defined as * = e*^ ''''/ {d—p)\. The classical 
gauge invariance of (|1.1| ) is given by 

6A = DAo , 

6B = -i[Ao,B] . (1.3) 

Using the field equations (which is equivalent in the path integral to integrate over B) 

F = , (1.4) 

DB = , 

we get the standard YM action 



Sym 
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\Jtt{FA*F) . (1.5) 



Note that the relation between the formulations (1J-) and ( |1.5|) may be also undertood as 



a duality map; however off shell B does not satisfy a Bianchi identity and this fact will 
be connected with the introduction in the theory of magnetic vortex lines. We remark 



that the short distance quantum behaviour of ( |1 . 1| ) and ( |1.5|) are the same as it has been 



exphcitely checked W2 . 



We remark that the action functional in (|1 . 1| ) defined by 




(1.6) 



is known as the 4D pure bosonic BF-theory and defines a topological quantum field theory 
0, |10[- The form of (|1 . 1|) was understood in |^ as an indication that the bosonic YM 
theory can be viewed as a perturbative expansion in the coupling g around the topological 



this perturbative picture around BF-theory we check in this paper the area law behaviour 
for the Wilson operator. 

The plan of the paper is the following: in section 2 we introduce an explicit analytic 
realization of the 't Hooft magnetic variable in terms of the A, B fields and check its 
properties at the operator level using the canonical quantization. In section 3 we compute 
its expectation value in the abelian projection gauge and verify the perimeter law. In 
section 4 we verify the area law for the Wilson loop operator. At last, in section 5 we 
draw our conclusions. 

2 A Color Magnetic Operator 

In this section we define a gauge invariant non local operator M(C) = M(E, C) associated 
with a fixed orientable surface S in (our base manifold) and with a suitable choice of 
a closed contour C on S, which gives an explicit realization of the 't Hooft loop. Part of 
this section was already investigated in by one of the present authors. 

The presence of the Lie-algebra valued two-form B field in ( p. . 1|) allows the definition 
of the observable gauge invariant operator 



where k is an arbitrary expansion parameter, x is a fixed point over the orientable surface 
S G M'^ (we do not integrate over x) and the relation between the assigned paths 7, 7' 



pure BF theory (|T]6|). This procedure has been called topological embedding |]Tl|. Using 




(2.1) 
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over S and the closed contour C is the following: let S to be a piecewise linear (PL) 
approximation of S by plaquettes (any two dimensional topological variety admits a PL 
decomposition). Then the closed path C is given by a succession of subpaths formed by 
the 1-skeleton of the 2-PL manifold S. In other words C starts from the fixed point x, 
connects a point y G S by the open path and then returns back to the neighborhood 
of X by 7^j;, (which is not restricted to coincide with the inverse {;^xy)~^ = lys)- From the 
neighborhood of x the path starts again to connect another point y' G S. Then it returns 
back to the neighborhood of x and so on until all points on S are connected. The path 
Cx = {7 U 7'} is generic and we do not require any particular ordering prescription as 
it is done in similar constructions devoted to obtain a non abelian Stokes theorem p3| . 
In figure 1 we show a typical path C connecting four points xl, x2, a;3, xA on a surface E 
given by a plane. 



Figure 1 

Of course the quantity (|2.1|) is path dependent and our strategy is to regard it as a 
loop variable once the surface S is given. In Eq. ( |2.1|) Hol|(7) denotes the usual holonomy 
along the open path 7 = '^^y with initial and final point x and y respectively, 

Hol|(7) =Pexp(i rrfx'^A^(x)) . (2.2) 

J X 

Given the finite local gauge transformations 

Hol| ^ g-\x)Y{o\lg{y) , (2.3) 
B{y) ^ g-\y)B{y)g{y) , (2.4) 

it is immediate to prove the gauge invariance of (|2.1| ). A comment on the geometrical 



meaning of the operator M is in order. If, for example, we take the surface E to be a 
torus we may define M(S, C) in terms of the parallel transport operator of loops |^ given 
by PexpjzA;/^^ B^{Ci)dx^^}, where B^{Ci) = §(^^{HolB^yHol')dx'' is a connection in the 
loop space and C is a linear combination of the two fundamental cycles Ci and C2 of the 
torus. 

By adding suitable bosonic vector and ghosts fields to the BF-YM theory, the "topo- 
logical" (we will clarify later what we use quotes here) symmetry of the pure BF action 
O, given by 

5A{y) = , (2.5) 
5B{y) = DA,{y)^d^y)A,{y)+t[A{y),A,{y)] , 



could be extended to the BF-YM theory. Then by using the well-known identity [jH 

d(,)(Hol|(7)Ai(2/)Hol^(7')) = Holi(7){rf(,)Ai(y) +^[A(i/),Ai(y)]}Hol^(7') (2.6) 



one may prove the invariance of ( p.l|) under the transformations (|2.5| ) up to boundary 
terms for non-compact S. However these boundary terms may be neglected in the renor- 
malized theory if mixed Dirichlet-Neumann boundary conditions for the quantum fields at 
the boundary 9(M^\S) = — 9S are used, in a way similar to the quantization of systems 
over a non complete space time (e.g. Casimir effect). 

Let us give a geometrical interpretation of the "topological" symmetry (|2.5| ). The 
holonomy dressed field Holi?Hol' can be interpreted as a colored charged field that trans- 
forms covariantly under the global colour symmetry g{x). Accordingly the transformation 
( p.5| ) represents a local change (with parameter HolAiHol') of the colour charge [|1^ . 



In the following the path dependence of ( p.l|) from the closed contour C will be crucial. 
Indeed, we shall see below that different choices of C for fixed surface S, we call it a 
framing (C, S), give different color magnetic fluxes in QCD. 

The classical BF-YM action ( [1.1|) is linear in time derivatives, and hence it is easily 
cast into its canonical form Let = x M and let t be a coordinate labelling the 
different A^^ surfaces. Let then be any vector satisfying ^^d^jt = 1. The appropriate 



notion of time derivative is then given by the Lie derivative 

L^A;: = dMin + Cd.Al - Cd.Al = D,{AX) + CF., . (2.7) 

one recovers the usual definition of time derivative. Time 
derivatives act only on the space components of A^. Therefore the conjugate momenta 
of y4g,-Boi are zero and these fields appear as Lagrangian multipliers since they may 
be eliminated through the classical constraint equations. In order to have first class 



constraints in the Hamiltonian formalism, the field content of the action (|1 . 1| ) must be 
enlarged in such a way that the resulting action enjoys the same topological symmetry 
( p.5|) of the action ( |1.6D . We thus introduce a Lie valued vector field rj which transforms 
under the symmetry ( p.5| ) as 61] = Ai and enters the action (|1 . 1| ) by substituting 
{B-Dr^f m- 

The 3+1-action becomes 

Sbf-ym = \ j dt j^, Sx{2ie^^^B1^Al + ^[i?o", - {D,r^^ - Ar/o")]' 

- ^e'''Bl,F^,, + iAy=\D,B,,r + ^^[B^^ - (Ar/J - D.r^t)?} (2-8) 

and the classical constraints read 

Ci[r(f)] = ^l^j'xiT^e'^\DiB,kr = , (2.9) 

C2HX)] ^ I f d'xvm^-^e'''F^l)=0 , (2.10) 

where a is the Lie algebra index, i,j, k = 1, 2, 3 and t"-, vf are arbitrary Lie algebra valued 
fields on A^'^. 11? is the canonical momentum conjugate to rj 

n- =y«-5or]f + z[Ao,r/,r-2[A„r/or} , (2.11) 

and 

xt = tU,kB';, (2.12) 

is the canonical momentum conjugate to ^4°. 

The appropriate gauge for the Hamiltonian formalism is the temporal gauge, which in 
this contest reads 

^0 = 0, (2.13) 



B^,, = . (2.14) 
fixes the conventional gauge symmetry and ( |2.14| ) fixes the topological symmetry. 



The Poisson brackets are defined by 

{A^{x), X%y)ho=yo = ^S''%,6^'\x - y), 

K(f),n5(y)} = -y), (2.15) 

and yield 

6At{x) = {Ci, A^ix)} + {C2, A^{x)} = -2iD,t\x) , 
SBt^ix) = {Ci, Bt^{x)} + {C2, Bt^{x)} = -2i[r{x), Bi,{x)Y 

~2i{D,v^{x) - Djv^ix)) , (2.16) 
6vt ={CuVm} + {C2,vm} = -^[rix),v^ix)r-vtix) . 

The constraint Ci generates the usual non abelian gauge transformations for A'^, rj^ and 
-Bj"-, and C2 generates the classical topological symmetry ( p.5|) . The canonical quantization 
of ( p.8|) is obtained by replacing the Poisson bracket {■, ■} in ( p.l5|) with the commutator 



■] and promoting the classical constraints Ci,C2 to operator valued ones 

Ci\A,B,r]> =0 , (2.17) 
C2\A,B,T]> =0 , 

where \A,B,ri > is called a physical state if it satisfies the above quantum constraints. 
We now show that M{C) = M{C, S = 5*^) generates a local singular (or equivalently a 
multivalued regular) gauge transformation, Qc{x), along C: this is precisely the defining 
property of the 't Hooft color magnetic variable 0. 

We start from the classical quantity associated to ( f^.lD and assume that S, and hence 
C, is small, so that one may consider a Taylor expansion of ( p.l|) around A; = 0. By 
using the Gauss theorem, the gauge ( p.l3| ), ( |2.14| ) and the identity ( p.6| ), which in our 
case becomes 

Tr[ J ^da'^iRoWijRoY)] = T^AJ ^ d^xe'^''di{}io\Bjk^o\')] = 

= Tr[/ d^x(Eo\ e'^''D,Bjk'Ro\')] = Tr[[ d^x}ior\C^)e'^''DiBjk\ , (2.18) 

JM3 JM^ 



one gets at lowest order in k 

M(C) ~ Tr{]l + ifc / d^x[Ror\C^)e'^''D,Bjk{x)]} , (2.19) 

where is a small 3-ball with boundary S'^ : dM^ = S"^ and Cx is defined as = 
7^ U y^. We now expand also Hol""*^ at this order 

Ror'iCx) - {H - lifci^)) , M^)={dy'A . (2.20) 



Substituting ( p.20|) in (|2.19|) and using the definition of the constraint Ci, Eq. (|2.9| ), we see 



that M{C) becomes the generator of the classical infinitesimal local gauge transformations 
with parameter kipc{x): 

M(C) ~Tr{Il + A;Ci[</.c(5^)]} , (2.21) 

where Ci = 2Tr[Ci] in the temporal gauge 0. 

To extend (|2.21| ) to the quantum level we substitute the Poisson brackets with the 



quantum mechanical commutators and rescale the geometrical fields to the physical ones, 

Any) = 9AT{y) , (2.22) 
[My), Bl{x)U=y, = j5'^%rs5^'\y - x) . (2.23) 

Owing to the cyclicity properties of the trace, ( p.l9| ) must be understood as symmetrized 
in DB and Hoi. Then, when (^.211) becomes operator valued, the ordering procedure 



required by the quantization and enforced by (p.l7|) , implies that we can substitute the 



product of fields with their canonical commutator p9|] . 

Therefore using the Gauss theorem of ( p^.lSD in ( [^.19| ), we express the divergence of 



the B field as the fiux through the surface S"^ and obtain 

M{C)\A >^ Tr{l + 2ik i dy' f da\^^eirs5^^\y - x)5"^T„Tfe}|A > . (2.24) 

^ Notice that M{C) is the non-abeUan generaUzation of the observable carrying non trivial magnetic 
charge studied in Ref.||l^, |l8j. Indeed the variable tt' defined in eq.(14) of Ref.|]l^ and in eq.(lO) of 
Ref. [T^ is the conjugate momentum of the vector potential which, in our first order formalism, is to be 
identified with the B two- form. Their monopole magnetic charge becomes, in our formalism, 4>c{x). See 
also next chapter. 
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Notice now that 



(f dy' (f dx''eiri^(—J-^. 
An Jc Jc oxi\\y — x\ 



d 



dxi \ \y — x\ 



Link((:7, C') 



(2.25) 



where S = S' U S", (9S' = C and C encloses the singularity of (|2.25| ). Here we have 
used the formulae 

d d 



i\\y-x\J 



dx^ dx_ 



dx^'x'' 



c 



(2.26) 



and the definition of the so called linking number between two closed loops, 
Link(C, C) 



/ dy' (f dx'eiri^i ^ 
Air Jc Jc oxi \y — x\ 



(2.27) 



irl ' 



{y - x)^ 



Att Jc jc \y — x\ 

where \y — x\ is the usual Euclidean distance between two points. Notice that in our 
construction the closed contour C is a framing contour of C, and hence the above linking 
number is the so-called self-linking number of C [pO|] , 



sLink(C) = Link(C, C'),^o , 
C^{x\t)} , 

C' = {y\t) = x'{t) + en'{t) ; e > 0, \n'\ = 1} , 



(2.28) 



where is a vector field orthogonal to C. The integral (|2.28|) is well defined and finite, 
takes integer values and equals the number of windings of C around C. 

Putting the above formulae in (|2.24|), we find 



M{C)\A> = \A + D^c>= 

~ Tr{Il(l + 2iA;c2(t)sLink((:7))}|A > , 



(2.29) 



where 



<5"XT6 = C2(t)Il 



(2.30) 



and 



C2{t) = (2.31) 



for the fundamental representation of SU{N). Eq. (|2.29| ) implies that M(C) generates 
an infinitesimal multivalued gauge transformation Lpci^)- Whenever x = {x*} traces a 
framing contour C of the closed curve C that winds n = sLink(C) times around C, M{C) 
creates a magnetic flux § 

$c = ^^^^^sLink(C) . (2.32) 

At this order of approximation the finite multivalued gauge transformation ^(^[x] 
generated by the action of M{C) over some state functional is given by 

M{C)\A{x) >= \Qc^[x]{A{x) + id^)nc[x] >~ Tr{e*3*cil}|A(f) > . (2.33) 

Because of the multi- valued nature of f2c[3^]) one has that 

Qc[xf] = Z^M^i] , (2-34) 

where Xj = x'(t = 1) = x'(t = 0) = xl is the base point of the loop C parametrized by 
C = {x\t) : < t < 1} and Z^^ is 



3(^9*0) ]x = gifcH^sLink{c)^ _ ^^2.35) 



Since A^^ = ^^^{A + d)ilc should always be single valued, Z^^ must be in the center 
of SU{N). To recover the standard form of the center, we normalize the free expansion 
parameter as k = 2tt/{N'^ — 1) and require a special framing for C, namely that n = 
sLink(C) G [0, . . . , — 1] . With these normalizations the form of the color magnetic flux 
is given by 

27rn 

*e = ^ . (2.36) 

Given the properties of M{C), its commutation relations with the Wilson line operator 
W{Ci) are easily deduced §, 0, |21|]: 

M{Ci)W{C2) = e''^''^^'^^''^^^''^^W{Ci)M{C2) . (2.37) 
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The previous equation is the well known 't Hooft algebra. The non trivial commutation 
relation between the operators M and W set their duality correspondence in the sense of 
the order and disorder operators in statistical mechanics. 



3 Computation of < M(E, C) > 

The purpose of this section is to compute the average of the BF-observable M(S, C) 
defined in Sec. 3. Namely we define the dual loop functional by the normalized connected 
expectation value 

< M(C) >,„„„^ ^^^^ , (3.1) 

where 

M(C)^M(S,C) = Tr[exp{27r^g^? |^ rf2yHol|(7)5(l/)HolJ(7')}] , (3.2) 

2 

<1> = j VAVBexp{- JtiIiB AF + A*B]} . (3.3) 

With respect to the notation employed in the previous section, we have defined the ex- 
pansion parameter k in units of the bare color charge g with a suitable normalization: 
k = 27rqg. Q 



Due to its SU (N) gauge invariance, (|3Tl| ) must be gauge fixed. Following Ref. 



we choose the abelian projection gauge. Before using this gauge we would like to make 
some comments. Magnetic monopoles are believed to play a major role in confinement. 
They appear as classical solutions of the Yang-Mills theory in the Georgi-Glashow model 
coupled with scalar fields transforming in the adjoint representation. Their classical mass 
is calculable, but turns out to be too big to allow for a magnetic Higgs mechanism. The 
Abelian projection gauge is supposed to generate monopoles of low or zero mass without 
introducing the scalar Higgs field. As monopoles transform as U{1) gauge fields, the 
idea is to isolate the U{1)^^^ degrees of freedom from those transforming in the coset 

SU{N)/U{1)^^^ with an appropriate partial gauge fixing leaving the U{1)^^^ maximal 

^ For the sake of generality we do not ask any flux quantization and, as a consequence, any fixed k 
value. After gauge fixing and a saddle point evaluation of the functional integral this quantization will 
naturally emerge. 
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abelian torus unbroken. The residual U{1)^~^ gauge invariance does not protect charged 
degrees of freedom from getting a massQ. Since confinement is a large scale phenomenon, 
in the crudest of approximations, these degrees of freedom can be discarded. Therefore 
the effective theory should contain only "photons" and monopoles. 

To implement this gauge we need a microscopic field or a composite of it, X, to 
transform in the adjoint representation of the gauge group SU (N) so that a gauge trans- 
formation can diagonalize it, 

X = g-'Xg = diagiXi,...,XN) . (3.4) 

The eigenvalues of such matrix can be naturally ordered Ai > A2 > ■ ■ ■ > A^v- When two 
eigenvalues coincide we have a line of singularities representing the monopole world-line 
By analogy with the interpolating gauges used in massive Yang-Mills |^2|, having 



chosen (for example) X = F12 and Yij = {Fi2)ij/{Xi — i 7^ j = 1,...,N, the 

proposed gauge condition is 

Y^'' + *A^'' = , (3.5) 

where the superscripts ch, stand for the off-diagonal and diagonal part of the matrix 
A = A'^ + A'^^ and is the covariant derivative with respect to the diagonal part of 
the gauge field. Interpolating gauges are such that for short distances the relevant gauge 
is * A'^^ = and the theory is renormalizable, while for large distances the gauge is 
ych ^ 0. For this reason the field Y needs to be adimensional and ^ has the dimension 
of the inverse of a mass squared. In standard YM theory, X can only be a composite 
of the F^^ (and its covariant derivatives) in order to transform according to the adjoint 
representation. But as the dependence of F^^, on the momentum is the same as the 
second term in the r.h.s. of (|3.5| ) the dominance argument is only possible if X is made 
adimensional. This makes this gauge difficult to apply to standard quantum Yang-Mills 
theory because Y is non-polynomial outside the tube around the monopole string and 
divergent inside it. 

^ This mass does not have to be physical and it can also be an infrared regulator. The decoupling of 
charged degrees of freedom is insensitive to the nature of the mass. 
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Quite remarkably, these problems are not present in the BF — YM theory due to the 
presence of the microscopic B field. The interpolating gauge can now be easily imple- 
mented choosing, for instance, X = B12: 

B"^ + D° * = . (3.6) 

Equivalently in our formalism we can diagonalize the two-form B on the surface S 

B =V-'BV = dtag{f]i,...,(3N) , Pi > (^2 > ■ ■ ■ > , (3.7) 
A =V-\A + d)V , F = V~^FV , 

and then use the background gauge condition D'^*A'^^ = in the renormalization program. 

According to the previous observations, in the computation of the expectation value of 
the observable M{C, S), employing the abelian projection gauge, in the large scale region 
we may neglect the massive off-diagonal degrees of freedom B'^^, A'^^. This approximation 
is often called Abehan dominance (@, @]). 

We would like now to interpret this approximation from a different point of view 
which will turn out to be very useful in the following. The existence of monopoles in the 
Abelian projection gauge is due to the compactness of the f/(l)^~^ group and it is related 
to the existence of non trivial topological objects for the entire SU (N) theory. So, in this 
contest, Abehan dominance can also be seen as the prescription of having reducible gauge 
connection which, as we will find later, will also be singular. In the simple SU{2) case 
this leads to a connection 

. 1 /a \ , , 

The reducibility of the SU{2) gauge connection, implies that the gauge bundle is split 
and thus requires the existence of a positive definite first Pontrjagin class and intersection 
number (which we will define later on). This fact, in its turn, implies the absence of anti- 
self-dual harmonic (closed) two forms |2^. Later on in this chapter, we will comment 
on the possible relation between this split bundle and the moduli space of instantonic 
classical solutions. 

In the general SU{N) case we now rewrite the functional integral in terms of the 
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variables 

ai=[A%, z = l,...,iV, ^«i = , (3.9) 

i 

and 

f^ = [F% = d[A% 2,j = l,...,iV . (3.10) 

It is then convenient to rescale the previous geometrical fields , to the physical 
ones A^ — > gA^, ' Furthermore, we replace the surface integral in Eq. ( p.2| ) 

with the integration over the so-called Poincare dual form, cj^, of (the homology class of 
) the surface S By definition, tus is closed. Moreover, choosing an orientation of the 
four manifold and remembering the absence of anti-self-dual harmonic two forms (which 
is imposed by requiring the existence of a gauge connection of the type of ( |3.8| )) cug is 
chosen to be self-dual i. e. *uy, = cus and with the property that (up to gauge one-forms) 
for a generic two-form t 



/ ujj^At . (3.11) 

In a local system of coordinates {x, y, u, v) on S, so that S is given by the equations 
X = y = 0, the dual form uj^ can be taken as cus — 5^'^\x, y)dx A dy and normalized as 



iV(S) 



5^'^\x,y)dx Ady =1 , (3.12) 



where iV(S) is the transversal tubular neighbourhood on the surface S. 

Let us now compute the average of the 't Hooft loop operator M{C) in the above 
abelian projection scheme. 

The functional integrations over the a^'s should be constrained by Yl,i = 0. In 
order to perform the calculations without imposing the previous constraint, we extend 
the SU (N) gauge group of the BF-YM theory to U (N) ||2^ . More precisely we identify 
the fields with the Cartan generators of U (N). Then we can recover the original SU{N) 
BF-YM theory by gauging the spurious U{1) symmetry, generated hj a^, associated to 
the standard embedding SU{N) U{N)/SU{N) = U{1). Finally the abelian gauge 
fields ai, ttAT-i associated to the Cartan subalgebra of SU{N) should be gauged to 
remove the remaining U{1) gauge degrees of freedom. 
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Using the Abelian dominance approximation, the 't Hooft loop operator M{C) be- 
comes 



M(C) =Tr{Oy(C)} , (3.13) 
Oij{C) = 6ij exp{i27Tq J ^ujj: A f3j[cos{g <j) aj) + isin{g<j) ctj)]} ■ 

As we have aheady observed, the gauge transformation in ( |3.7| ) is singular 0, and the 
singularities occur if two consecutive eigenvalues of B coincide. If /3j = we shall label 
such a point by x*-*-*. Moreover is a two-form and hence its support is a two-cycle in 
the four manifold. Thus, in order to avoid the previous singularities x^^\ we must remove 
a ball S'^{x^^^) of radius e — ^ 0+ at each point x^^^ of the base manifold. As a consequence 
of that, we get a non trivial Stokes theorem {gai is the geometrical U{1) connection ): 



g f at = g / dai = g dai = 271- goi (3.14) 

Jc Js : dS=C Js^ix(i)) 

{S = S'US'^{x^^^)), since the two-cycle S is not homotopically trivial (it contains S'^{x^^^)). 
In ( p.l4|) the magnetic charge is given by the first Chern class, ci(L), of the line bundle L 



with fiber S'^{x^^^). Since we have in principle (A^ — 1) connections U{1) whose curvatures 
are supported near the points x^*^ of S (x^^-* is identified with x*-^-*) and the total first 
Chern class is an integer (ci 

= (AT - l)ci(L) e Z), we obtain that gqi satisfies 

the Dirac quantization condition rule 

3^1^ = ' ^^ = ±1, ±2, - , V^G[l,Ar] . (3.15) 

In other words the appropriate homotopy group is 'n.2[SU {N) / U [1)^ = Z^~^, as ex- 
pected. 

In the following we shall show that the effective field equations associated to the strong 
coupling limit of T{C) = — In < M{C) >conn admit as critical points magnetic monopole 
strings - which we shall construct explicity - with winding numbers nj G Z at the singular 
points x*^*\ The key point will be the identification of the arbitrary expansion parameter 
q in ( p.2|) with the magnetic charges g^, i.e. we will set 

qi (X q (X - . (3.16) 
9 
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We start now the computation of the magnetic order parameter ( |3.1|) in the g — 0, which 
will be identified by as the strong coupling limit (7 — > oo. In this limit the operator 

( p.l3| ) can be approximated by | 



M(C) = 5]0,,(C)~iVexp|^y'a;sA^Acos(^j^«*)} • (3-17) 



Eq. (|3.17| ) is obtained by taking into account that, according to our previous observations, 



there are two types of monopole configurations for each i, depending on the sign of the 
magnetic charges g^. Indeed we have that {qi} = {g|} = {ql' > 0,q~ < 0}. Hence it is 
resonable to define the sum over the field configurations as J2i = J2i,s^ where i E [1,N] and 
s = ±. Notice that < M(C) >conn in the strong coupling limit is the partition function 
of N compact QED's or 4D Villain models (e.g. see ||^). 

With all these approximations taken into account the form of the magnetic order 
parameter in the strong coupling region becomes 

N r I ^ r 1 

< M{C) >c^n-^Y^ j VaiVpi exp{-- ^ J A *A + i(3, A [dai 

+ ^^scos(^^aj)])} = J Vaiexp{-^Y1 J [d^i + ^^^cos{g (^^ai)]'^.18) 



where the square of a form t means t A *t. We now split the gauge fields as ai = ai + hQi, 
where the quantum fiuctuations Qi must be gauged (e.g. by a covariant gauge condition) 
and the are singular classical configurations. Postponing for a while the discussion 
of quantum fiuctuations, we concentrate on the semi-classical contribution to the path 
integral which is 

/ cje A cjs' cos{g 

ic ' Jc 

In the derivation of the above equation, we have exploited the following facts: 



< M(C) AT exp{-^^^ j A^s'Cos(c/£ai)cos(^£^a.)} • (3.19) 



partial integration is not allowed on the ai due to their singular behaviour [25], 



^We note that the appearence of the effective expansion parameter q/N in ( ^.17 ) suggests that the 



same conclusions of this section can also be obtained in the large TV-limit (holding q fixed) in which the 
SU{N) and the U{N) gauge theory are the same Q. Thus, in the q —^ limit, we can consistently drop 
the constraint ai = 0. 
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• ttj obeys the (monopole) equation 

4:7rq f 
* dtti = -r^cjs cos(5f f ai) , (3.20) 
A* Jc 

• the absence of electric currents in the model, 

• the self-duality and closedness of cus. The above properties imply the absence of 
terms linear in Qi. 

Equations of the type ( p.20D , appeared already in the study of the duality properties 
of gauge theories and 4D manifold invariants P7| . 

The r.h.s. of ( p.20|) can be written in terms of spinor fields for a direct comparison with 



Ref.|p7||. If M'^ is a spin manifold any self dual two- form B~^, satisfying the conditions (in 
spinorial notation) B^^^ A 5^^^ = and B^^^ /\B\q ^ 0, can be represented in the form 
[B^)"^^ = ^{MaMb + MbMa), where M is a two component spinor field. Furthermore 
^^uV'ilu:^'^AB = ^abj "where t] are the 't Hooft symbols and a the Pauli matrices. In 
terms of the space-time indices 5^^, = MF^^M, with F^,^ = | [7^,71/]. Using the equation 
of motion for the B field, M is seen to obey the Dirac equation relative to the case 
of manifolds which are Kahler and spin as discussed in section four of Ref. ||27||. This 



connection with Ref. |2^ , in which the infrared limit of the supersymmetric N = 2 theory 



is described in terms of abelian gauge fields and monopoles, should not be surprising. In 
fact, if we take the attitude (that will be developed in the next section) that the term 
g'^/4:B A *B is a perturbation of the "pure BF" theory (|1.6|) , this perturbation theory will 
be performed around the vacuum of the theory which coincides with that of ( |1.6| ). This 
requires to gauge fix the local symmetries ( p.5|) , (|1.3| ) of the action ( |1.6|) by introducing 



appropriate ghosts |10 



We would now like to explain in which respect the pure BF theory (with action (|1.6| )) 
deserves to be called "topological". In first place the action ( |1.6|) is metric independent. 
Moreover, after some work, a BRST charge , s (nihilpotent on-shell) and a vector charge, 
6^, operator can be defined for the quantized theory [jlO|. The commutator {s,6^} closes 
on translations on-shell and this shows that the stress energy tensor of the theory can 
be written as a BRST commutator. Introducing external sources coupled to the BRST 
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variations of the microscopic fields a new linearized BRST charge, Q (nihilpotent off-shell), 
can be defined. Thus the BF theory is proved to be "cohomological" [PUI with respect to 



the f2 operator but it is not clearly equivalent to the topological theory defined in [2S] 
since the BRST charges of the two theories are different. 

Independently from these considerations, the symmetries (|2.5| ) must be anyway gauge 
fixed, which can be done using a covariant gauge for the field B. 

We now look for an explicit solution of (|3.2CI|) under the constraint cos{g§^di) = 1 
which is consistent with the Dirac quantization rule ( |3.15| ) if we set = m{N — 1), m G 
^5 1i = Jil- We then get 

mN 

^9 = ^ ■ (3.21) 

Given an arbitrary and topologically trivial surface S the dual Poincare form can be 
locally written as 

= J^da,MS^'\x-y) , (3.22) 

where are the coordinates on the surface, are the normal vectors of the bidimensional 
surface S and a = 1,2. Putting ( p. 221) in (|3.20|) we find the classical background solution 



31, 32 



al{y) = % d'^P^^^^f.e'^Xt'^d^j-^-^ , (3.23) 

where P'^ = (1 + *)/2 and, as before, x*^*-* is the location of the singularity at the in- 
tersection of the surface S with the i-th monopole world-line and S locally appears as 
the "Dirac sheet" ^j. The curve C of (|3.20|) clearly winds around x*^*-*. Eq. (|3.23|) is 
a magnetic vortex-line or Abrikosov-Nielsen-Olesen (ANO) string for each i G [l,iV] Jill) 
obtained here without the use of the Higgs fields, but introducing the disorder effect of 
the M{C) operator. [] We stress that these configurations do not correspond to 't Hooft- 

Polyakov like monopoles but are rather singular Dirac monopoles [Q. The bidimensional 

^ In Ref . p^, 1^ the same solution is found in the abehan Higgs model with aj^ — qidf^Xi ^^'^ the Higgs 
field is parametrized as ip — \ip\ exp(ix). 

In Ref.js^ these kind of configurations are shown to describe massless monopoles in 4D compact QED, 
giving rise to a planar model like thermodynamics. 
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reduction of the solution ( |3.23| ) (obtained fixing the surface E) coincides with eq.(31) of 
Ref.p^l (near the singularity) which is a solution of the dimensionally reduced monopole 



equations of Ref . [ 

Finally, the classical contribution is 

So = — I^uj: Aus^ cos{g j ai) cos{g ia'i) = (5(S, S') . (3.24) 

For ( p.24| ) to make sense, the two cosines must be different from zero which implies that 
the paths C, C must encircle monopole flux tubes. There are now two options: 

• 9§c0^ii^ 9 §c' 

• 9§c^i = 9 §c' ■ 



In the first case, Q is the intersection number |35 



Q(S,S')=^^^ws[c]Aa;s'[c'] , (3.25) 

which is an integer (it is a four dimensional topological invariant) characterizing the self- 
intersection of S with itself. This number is given once the surface S and the ambient 
space are assigned. 

In the second case the closed curve C = {y^{t)} is a framing contour of the closed 
curve C = {x^(t)}, i.e. iff it happens that 

y^(t) = x^(t) + en^(t) , e ^ , \n''\ = 1 , (3.26) 

where n^{t) is a vector field orthogonal to C. In this case Q(S, S') becomes the self-linking 
number of C given by ( |2.28| ). 

Then we get that 

Q(S, S') = sLink(C) . (3.27) 



From a physical point of view, we may define sLink(C) in ( |3.27 ) as 



sLink(C) = , (3.28) 
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where p plays the role of unit of lenght and L{C) is the perimeter of the loop C. We may 
understand ( ^.281 ) by a n- vertex polygon discretization C ^ CpL of the loop C. Indeed by 
the very definition of the self-linking number one gets that sLink(CpL) = n = L{Cpl)/1, 
where / is the lattice spacing. ( p.28| ) is obtained in the continuum limit n ^ oo, / — 0, 
since in this limit L{Cpl)/1 ~ L{C)/p. Thus, in this picture, the perimeter law comes 
out from the arbitrary framing dependence of the path C in the color magnetic operator 
M{C, S). In turn, the framing dependence can be seen as a point splitting regularization 
due to the non-local nature of the operator. In our construction sLink(C) in 4D cannot 
be slinked - in general the linking between two arbitrary curves in 4D has no geometrical 
meaning - since C and C belong to the same surface. 

Let us now discuss quantum fluctuations. If the effective theory for large length scales 
is a U{1) type theory, for short length scales the charged degrees of freedom cannot 
be discarded anymore. Let A be the length scale separating these two regimes and let 
us divide the gauge field according to a background field prescription: A"- = A"- + 
where a is a gauge index. Moreover the gauge field can be written as the sum of two 
parts giving the contribution of the large length scale sector (> A) and the short length 
scale sector (< A); the functional integration over the gauge field factorizes in a way 
compatible with this separation. Moreover let, for the sake of simplicity, the gauge group 
be SU{2). For length scales bigger than A we take A^ = a + which is the usual U{1) 
prescription. For length scales smaller than A we take A"" = S"'^a + Q"", i.e. we continue 
the classical solution into the small scales region where the quantum fluctuations coming 
from the charged degrees of freedom cannot be discarded. The expectation is that the 
small scales behaviour is insensitive to the classical solution according to the background 
field method. Performing the functional integration over the quantum fluctuations leads 
to a double contribution, in complete analogy with the saddle point evaluation around an 
instanton background PB|. 



The first contribution is given by a ratio of determinants given by 

(3.29) 



Det'{-Lo) 
Det'{-L) 



Det{-d^) 
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where D = d — iga, D'^ = D'^Df^,, 

L = DH^,-{l-^)D^D, + 2igFl^{x) , (3.30) 

and Lq is given from L evaluated around the trivial background. ^ is a gauge pa- 
rameter usually chosen to be one. In (|3.29| ) the determinants are primed to remind 



the reader that the contribution of zero modes is omitted and that the determinants 
are regularized. 

• The second contribution is given by the Pauli-Villars regularization of the determi- 
nants and it amounts to a scale /i (which is the Pauli-Villars mass) raised to a certain 
power which is given by the dimension of the moduli space of the classical solution. 

Let us now proceed with the evaluation of these two contributions. Using the self- 
duality property of our classical solution, the ratio of determinants ( |3.29D can be written 
Det{—d'^)/Det{—D'^) |3^. This ratio has been evaluated in Ref . |p7|| using the 



as R 



heat kernel method in the case of an SU (N) gauge group but it is easy to generalize this 
result to our case too. We now give a brief outline of this computation to justify the 
previous statement. 

To be mathematically well-defined, the split connection (p.8|) must be considered either 
to be regular in 3?^ or to be singular in S"^. In our case we found explicitely a singular U{1) 
connection and we can go from one point of view to the other by removing the singularity 
locus spanned by the a;*^*^ in (|3.23|) . Indeed, by removing a disk from S'^ we in fact obtain 

The computation we want now to perform is intended to probe the region around 
the singularity by evaluating the operator R obtained by expanding around the classical 
solution given in (|3.2(]|) . This computation is carried on a space-time bounded by the 
length scale A which can be taken to be the compact manifold 5*^ after a conformal 
transformation, to allow a proper definition of the eigenvalue problem and to avoid infrared 
problems. 

The starting point of the computation is the introduction of the Riemann zeta function 
({s) built out of the eigenvalues of the operator —D^. The logarithm of the ratio of the 
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determinants of the self-adjoint operators R is given by the variation, under conformal 
transformations, of the derivative of the Riemann zeta function C'(0) continued to zero 
value of its argument. In turn this function is given by the residue of 

2 / t'-^i:i[G{x,x,t)5oj]dt (3.31) 
Jo 



at s = [0. In ( |3.31| ), 6uj is an infinitesimal conformal transformation and 

1 (^_„)2 OO 

G{x, 2/, t) = ^^^e-T Yl ^n{x, y)r (3.32) 
satisfies the heat equation. 

Making use of (|3.32|) , (|3.31|) can be cast in the form 



R = ('(0) = At / Tra2(x, x)Sujd^x . (3.33) 
Svr"' J 

Plugging the expansion ( p.32| ) in the heat-kernel differential equation we obtain a set of 
recursion relations from which the values of the coefficients a„ are easily extracted. Since 
the only manipulation needed in solving such recursion relations is the commutation of 
covariant derivatives, the final result is easily generalized to any gauge group 

(3.34) 

The factor 1/4 comes from the normalization of the gauge group generators according to 



The contribution coming from the regularization of the zero modes is obtained once 
the dimension of the moduli space is computed, according to Ref.E^, to be 



dimM = ci{Lf = ci(L) A ci(L) = / {da)^i,{x){da)^i,{x)d'^x . (3.35) 



Putting together the classical result ( |3.24| ) with the quantum fiuctuations, we find that 
the bare coupling g can be substituted by its renormalized expression and that (|3.19|) can 
be written as 

S^^^(i?-(-f)-H^^(5;^-..MM)) , (3.3, 

where /3i = y is the first coefficient of the SU{2) beta function of the non-abelian Yang- 
Mills theory. 
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4 The Average of the Wilson Loop 



In this section we shall compute the average of the Wilson loop and find an area law 
behaviour for its leading part. Furthermore, in our formalism, the area law gets a nice 
geometrical interpretation: it is the response of the true QCD vacuum to arbitrary defor- 
mations of the quark loop C. 

The starting point here is given by the Wilson loop operator written in terms of the 
non abelian Stokes theorem (see e.g. ||13[ ): 



WniC) = Wn{^,C) = TrR{Pj,exp[i [ Hol?(7)F(a;)Hol^(7')]} , (4.1) 

where C = C = {'y{x) U7'(x)} was defined at the beginning of section 2 and P^, means 
surface path ordering. Wr is calculated with respect to some irreducible representation 
R of SU{N). The loop average in the fundamental representation R = t of SU{N) is 



where Sbf-ym was defined in (|1.1|) . 

Expanding in series the Wilson loop ( |4.1| ) we get 

< W^t(S,C) >= E < ^Trt^E / . . . / (2Hol(7)F(x)Hol-^(7'))" > • (4.3) 



We then use the identity 



-i!*Bt..F%FUx)=Ai ^--(e-i/*^-^M^) . (4.4) 



Performing a partial integration with respect to the functional derivative in (|4.4| ) we can 
replace, in the path integral, 

mol(7)e-i/*^°-^°F(3:)Hol-^(y) ^4Hol(7)e-i/*^°-^°( , 1, j Hor^(y) . (4.5) 

* B[x) 

The functional derivative acts only to its right on the exponential of the mass term 
—(7^/16 / B'^^B'^^, since Hoi does not contain the B field. We now need the identity 

)p<7 



V( - )e-i/^M.^2. = V(-— * )e-^/^M^^M. (4 6) 
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where V is the functional defined by 

6 . ^ r r ,. b 



^(?7B^)-^=4-l(^Hol(-,)(j^)Hor.(y)r . (4.7) 
Resumming the exponential series for the Wilson loop we finally get the "duality" relation 

<Wt{C) >,onn= — , (4.8) 

where 

2 

M;(S,C,C = 9S)=rr,[Psexp{-|- 1^ Hol?(7) * 5(x)Hol^(y)}] • (4.9) 

M^*(S, C) is the dual (in the sense that B — > *5) of the observable Mf(S, C) defined in 
(13) with k set to A; = ig^ jl. 

To calculate ( [4.8|) we expand perturbatively in g both the exponential and the holonomies 
which appear in the exponent of ( ^.9| ) p]|]. The first relevant contraction encountered at 
lower level is given in terms of < A * 5 >, which can be computed starting from the off- 
diagonal propagator < AB >. The latter propagator is the same for the BF-YM theory 



1^ and the pure BF (|1.6|) . Therefore we find 



<M;(C) e-^^^W^c/.<^*^>A(S) , (4.10) 



where C2(t) is given in ( p.31|) and A(S) depends on higher order integrations over E. 



In the limit g ^ ( p..!] ) reduces to ( |1.6| ) which is a topological theory. Making the 
perturbative expansion in g around the topological theory, A(S) in ( [4. 10|) is found to be a 



diffeomorphism invariant quantity, i.e. A(S) depends only on the knotting properties of S 
and does not depend on the linking properties of S (at least when our ambient space-time 
is S^) . If the dynamical surface S created by the quantum fluctuations of the true QCD 
vacuum is the unknotted surface Eq, one should flnd that 

S = So^A(So) = l . (4.11) 

The explicit calculation of A(S) is an open problem, but for the purpose of showing the 
area law behaviour its knowledge is not essential. 
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Consider now the integral in the exponent of ( [4.10| ), 

/ / <A*B>=< (f A f daf"'{x){*B{x)).^ > (4.12) 

Jc JT. Jc JT, 

=< j^A lj*dar{x)B{x)^, >= £dy^lj*danx) < A{y)xB{x))^, > , 

where *da{x) is the infinitesimal surface element of the plane S* dual to S at the point 
X G S. We may rewrite 

/ dy\*danx) < A{y)xB{x))^, >= I dy^ l^^"< AxB^. > , (4.13) 

recalling that cus is locally given by a bump function with support on S. Eq. ( [4.13|) is 



by definition the linking number between the curve C and the dual plan S* in x to S. 
Indeed the linking Link(C, S), whith arbitrary C and S, is defined by 

Link(C,S) = 1-1 dx^ I daP\y)e^p,SY^^ = A f dx^ [ da^^y) < A^(x)i?^^(y) > 

(4.14) 

In our case, by construction, Link(C, S*) ^ 0. The residual integration over S in (|4.13|) 
spans all the dual S* to S, yielding a contribution proportional to the area of S, 

<m(C)>.onn ~ exp{V^^^^^^/^Lmk(C,S:)} , 
/ Link(C,S;) ~ A(S) . (4.15) 

JT, 

We may get a better understanding of ( [4. 15|) considering a lattice regularization of S, i.e. 



S TjpL. In this case C runs over the links of the lattice, while S* corresponds to an 
element of the dual lattice through the plaquette centered at x. J-^^^ Link(C, S*) is an 
integer which counts the number of vertices of the dual lattice on Sp^, or equivalently 
the number of plaquettes Np of Sp^,. We may then write 

/ UnkiC,K)=Np = ^^^P^ , (4.16) 

where A{Tjpl) is the minimal area bounded by the quark loop C in the fundamental 
representation and a is the lattice spacing. Passing to the continuum limit a — >^ 0, 
Np = A(SpL)/a^ becomes A(S)//^ with / a typical scale in QCD which may be choosen 
as Aqcd- Therefore one may rewrite ( [4.15|) as 

< Wt{C) >conn-^ exp[-a{AQCD)A{E)] , (4.17) 
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where (t(Aqcd) is the string tension defined by 

- 1 

a{^QCD)^9l{-^)^%CD , (4.18) 

with Aqcd = 1/^- Here we have replaced, as required by the renormahzation of the 
theory, the bare couphng constant g with the running couphng constant at some QCD 
energy scale Aqcd- Of course we may change the physical scale, but it is well known that 
this corresponds to a different choice of the renormahzation scheme, the difference being 
only a finite renormahzation, leading to the same physical results. 

Existing data corresponding to energy scales between 4 and 100{Gevy can be fitted 
with Aqcd ~ 0.5 Gev, corresponding to [gj^{—p'^,AQCD)/4:Tr] between 0.2 and 0.4 for 
Nf = 4. Inserting these data in (|4.18| ) we find that a rough estimate of the physical value 
of the string tension can be given by 

Va{AQCD ~ 0.5Gev) ~ 316 - 440Mef , (4.19) 

which is of the same order of magnitude of the experimental value. 

Before concluding this chapter let us comment on the result we have found. It might 
seem puzzling that the area law behaviour comes out of a perturbative calculation. In 
reality our perturbative expansion is an expansion around the gauge fixed topological BF 
theory whose vacuum contains instantons. It is thus resonable to imagine that the non 
perturbative information comes out of the non-trivial vacuum structure of the quantum 
BF theory. Furthermore we remark that the computations of this section could also have 
been done along the lines of the chapter three, given the similarity between M(S, C) and 
M*(S, C). As a consequence the invariant S') in (|3.27| ) is replaced by Link(S*, C), in 
agreement with ( [4.15| ) and abelian dominance translates into ( [4.11D , i.e. lack of interesting 
knotting properties of S. 

5 Conclusions 

In this paper we have shown that in the framework of the BF formulation of Yang-Mills 
theory, 't Hooft views on confinement can be explicitly realized in the continuum formu- 
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lation. This has been possible because a non-abehan color magnetic operator has been 
identified quite naturally in terms of a non local observable (|2.1|) constructed out of the 



microscopic fields A, B. Then we have computed the average of this operator in terms of 
a saddle point approximation having used the abelian projection gauge. The implementa- 
tion of this gauge has consisted in performing the computation in terms of the large scale 
effective lagrangian obtained neglecting massive charged degrees of freedom. This point 
clearly awaits further clarification and an estimate of the contribution of the discarded 
sector remains to be evaluated. However, we do not believe that the perimeter law be- 
haviour in (|3.34|) , coming from the evaluation of the classical action, can be disrupted by 



quantum corrections. 

In any case the ratio of two operators of the type ( |2.1| ) with different paths C, C still 
exhibits the perimeter law since any possible divergence comes from shrinking to zero the 
size of the magnetic flux tube (ANO string). 

In the framework we have investigated, the perimeter law stems from the breaking 
of the BF topological theory (|r6D due to the g'^i:i{B A *5)/4 term in (|L1]). In [|g the 



expectation value of M{C, S) was computed in the pure BF theory looking for two-knot 
invariants i.e. for four dimensional topological invariants of knotting of surfaces. In that 
computation a classical action piece like (|3.24 ) did not appear because of the missing 



Tr[B A *B) term (which allows for gaussian integration) and only functions of the four 
dimensional topological invariant ( p.25| ) were obtained. It is thus reasonable that in the 
BF Yang-Mills appears a homotopy invariant of the curve C as sLink(C) which represent 
the effect of the winding of C around the monopole flux tube (which we can think closes 
at inflnity). The curve C in (|2.1| ) is connected to the holonomy of the gauge connection 
and contains the disks S'f (x(j)) C S. These disks become the observable Dirac surfaces of 
and are responsible for the perimeter law. An argument similar has been developed 
in for a BF Yang-Mills theory in two dimensions where the two form B becomes a 
scalar. 



Let us now remark on some properties of ( p.34|) which agree with common wisdom: 



the appearence of qr has been conflrmed in lattice simulations of SU{N) QCD The 
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"photons" Qi interacting with a plasma of condensed monopoles reproduce the same 
SU{N) running couphng constant required by the renormahzation group. From ( |3.24| ) 
we have a perimeter law behaviour for m 7^ and finite N as expected for the confining 
phase of QCD in the absence of massless particles. 

As far as the average of the Wilson loop is concerned, ( [4.15|) is the basic result. Indeed 
it implies that the area law for the Wilson loop average is a consequence of the response of 
the QCD vacuum to arbitrary path deformations of the loop C. This property is evident 
in our formalism, where the loop C = 9S may be understood as the deformation of the 
loop C. The existence of Link(S*,C) means that the expectation value < Wr{C) > is 
not invariant, as reasonably expected, under the change of the path-prescription C in the 
non-abelian Stokes theorem given by (|?T). In other words, ( |4.15| ) is a non-perturbative 
Ward identity associated to the path deformation of the static quark loop C. It is well- 
known that while the fundamental representation of the gauge group leads to confinement. 



the adjoint representation does not. This does not show up in ( |4.15|) due to presence of 
the coupling g which should be absorbed by the presence of a mass gap (from monopole 
condensation) in the propagator for the A and B field. 

Finally, observe that the results obtained in this paper are due to the non-abelian 
nature of the model. If we had started with a pure abelian BF-YM theory, the operator 
M would have been defined without the Hoi terms and, therefore, without the possibility 
to have a non-trivial monopole equation (|3.2CI|) . In the abelian case an additional Higgs 
field is needed in order to have a non-trivial magnetic variable M. 
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